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Abstract 

Based on the matrix Green's function method combined with hybrid tight-binding / density func- 
tional theory, we calculate the conductances of a series of gold-dithiol molecule-gold junctions including 
benzenedithiol (BDT), benzenedimethanethiol (BDMT), hexanedithiol (HDT), octanedithiol (ODT) and 
decanedithiol (DDT). An atomically-contacted extended molecule model is used in our calculation. As 
an important procedure, we determine the position of the Fermi level by the energy reference according 
to the results from ultraviolet photoelectron spectroscopy (UPS) experiments. After considering the 
experimental uncertainty in UPS measurement, the calculated results of molecular conductances near 
the Fermi level qualitatively agree with the experimental values measured by Tao et. al. [Science 301, 
1221 (2003); J. Am. Chem. Soc. 125, 16164 (2003); Nano. Lett. 4, 267 (2004).] 



1 Introduction 

Molecular devices are now playing an important role in nanoelectronic researches. The construction, mea- 
surement and understanding of the transport properties in molecular electronic circuits have drawn more and 
more attention. At present, a lot of theoretical work is concentrating on the sandwichtype structure of an in- 
dividual molecule between metallic electrodes [1-4]. Although current simulation methods can qualitatively 
describe some electronic properties of such simple systems, there are still a lot of issues waiting to be made 
understandable for the purpose of predicting, designing and controlling new molecular devices. For example, 
the methods used to calculate the coupling strength of electrodes-molecule interaction remain unconvincing. 
It's difficult to model the metal-molecule-metal junction accurately because of the uncontrollable interface 
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structures between molecules and leads. In addition, the position of the Fermi energy is another focal point. 
The common simple viewpoint is that the Fermi level of electrodes lies within the gap of molecular orbitals 
between HOMO (highest occupied molecular orbital) and LUMO (lowest unoccupied molecular orbital) [5]. 
Further calculation and experimental results are essentially needed to clarify the doubt of precise location 
of the Fermi energy. On the other hand, the importance of the modification of molecules due to the non- 
equilibrium situations and the determination on the spatial profile of the electrochemical potential across 
the interface also need more discussions. 

For device under equilibrium, the above problems reduce to calculating electronic structures of the 
molecules sandwiched between two metal electrodes, which is the first step to understand the transport be- 
haviors of molecular devices under any circumstances. Recently, the conductances of a series of gold-dithiol 
molecule-gold systems were measured near low bias by the tip of scan tunneling microscope (STM) [6-8]. 
Though the molecular junctions based on Au-S bond have been investigated extensively [9-14], our calcu- 
lation and further discussions according to the experimental results still can show more enlightenments for 
understanding the interaction between molecules and metals. Ultraviolet photoelectron spectroscopy (UPS) 
was referenced for the purpose of determining "line up" about the Fermi level [15-19]. UPS experiments 
can provide direct spectroscopic information about the gap between HOMO and Ep, which is one of the key 
points in the transport calculation of molecule-metal sandwichtype structures. In this paper, we present a 
first-principles calculation, which is based on the methods of Xue-Datta-Ratner [2,3,9], on the conductance 
of molecular wires including benzenedithiol (BDT), benzencdimcthanethiol (BDMT), hexanedithiol (HDT), 
octanedithiol (ODT) and decanedithiol (DDT). An atomically-terminated extended molecule model is used 
and our calculations arc in qualitative agreement with the experimental results [6-8] . We also try to analyze 
the influence from the derealization of molecular orbitals on the transport properties of molecular junctions. 

2 Methods 

We use the matrix Green's function method combined with hybrid tight-binding / density functional theory 
(DFT) for modeling molecular devices, which was firstly introduced by Datta et. al. [20,21]. Since what 
we now consider is the transport properties under equilibrium, the most important thing is to calculate 
electronic structures of metal-molecule-metal junctions. This problem is the generalization of the familiar 
chemisorption where a single isolated molecule is considered to be adsorbed onto one metal surface. The 
molecular device can be divided into two different parts: "extended molecule" region (the molecule and 
the neighbour metallic atoms on the electrodes perturbed by the adsorption of the molecule) and "contact 
region" (the unperturbed part of the electrodes) [21,22]. This separation provides a convenience to treat the 
two parts by well-established techniques independently The solution of the total device system is based on 
Green's Function method, which supplements the coupling of the molecule and electrodes via the concept of 
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self-energy. The retarded Green's function G R of the extended molecule is given by 

G R = (ES - F - Sf - S^)- 1 (1) 

where E is energy, S is the overlap matrix, F is the Fock matrix of the isolated extended molecule. F and S 
can be obtained through DFT calculation, which is performed using GAUSSIAN 03 program package [24]. 
The self-energy matrices Sf and Y> R can be written as 

Sf = TiVn and Z R = 4g R T 2 (2) 

where T1.2 is the coupling between the extended molecule and electrodes, and gi.2 is the surface Green's 
function of the electrodes. 

The computation of the surface Green's function of gold electrodes and the coupling matrices is described 
using tight-binding scheme with the parameters given by Papaconstantopoulos [23]. Here we suppose that 
the gold electrode is semi- infinite Au(lll) surface and calculate the Green's function layer by layer with 
an iteration procedure [21]. Then we can obtain the transmission function through the molecule by the 
following equation 

T(E) = Tr(T 1 G R T 2 G A ) (3) 

where 

r 1>2 = i(E* 2 - S^ 2 ) (4) 
Since what we concern is the conductance under equilibrium condition, the relation can be given by 

G= 2 -^T{E S ) (5) 

The density of state of the extended molecule D{E) can also be calculated by Green's function 

D{E) = — —Tr[Im(G R S)] (6) 



3 Results and discussions 



3.1 Optimization and calculation models 

In the following DFT calculations we have used the Becke's three-parameter hybrid functional using the 
Lee, Yang, and Parr correlation functional [25,26] (B3LYP). The basis set chosen for C, H, S atoms is the 
triply split valence with polarization functions, 6-31 1G** [27,28]. When the Au atoms are included in the 
calculations, two Los Alamos National Laboratory basis sets [29,30] are used with effective core potentials, the 
LanL2DZ basis set for the geometry optimization and the LanLIMB basis set for the transport calculations. 
The use of the pseudopotentials with relativistic corrections has been widely demonstrated to be a good 
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compromise with the alternative use of full-electron procedures, especially when we deal with relatively large 
numbers of heavy atoms. It requires less computational effort without loss of accuracy [31]. 

Under equilibrium condition, we have considered the electronic transport properties in molecular junc- 
tions consisting of benzencdithiol (BDT), benzenedimethancthiol (BDMT), hexanedithiol (HDT), octancdithiol 
(ODT) and decanedithiol (DDT). As the first step in our calculation, we optimize the molecules by a gold 
trimer-molecule-gold trimer structure, which is shown in Figure^ The optimization scheme is similar to our 
recent work on calculating 4,4-bipyridine connnected to gold electrodes [32]. The Au-S bond lengths of five 
optimized metal-molecule-metal structures are 2.32 ~ 2.33A All the optimized bond lengths are consistent 
with other theoretical data [33]. On the other hand, the calculation of Kruger et. al. also indicates that 
the gold trimer- molecule structure is stable after leaving from the gold surface [34], which supports our 
optimization scheme. 

Assuming the molecular structure remains unchanged during the pulling process by STM tip, we can 
use the above optimized structures in the following transport calculations. The contact geometry is always a 
debating problem in the construction of molecular junctions. Some previous work have already done careful 
discussions on this issue [34,35]. Because different adsorption geometries may correspond to dissimilar 
local minimal energy points [36], we can not determine the accurate contact by only a simple junction 
model. The work of Emberly's group shows that the conductance of gold-BDT-gold junction changes little 
when the geometries of contact vary near zero bias [35]. Their results indicate that an atomically-contact 
extended molecule model could be one of the reasonable qualitative descriptions for molecular junctions under 
equilibrium. Thus we model the junctions by inserting a single gold atom between the end sulfur atom and 
the gold surface at each contact. The inserted gold atom is settled above the center of 3-fold hollow site on 
Au(lll) surface and maintains the Au-S bond length calculated in the above optimized structures. Three 
nearest-neighbor atoms on the surface are included in the extended molecule, which is showed as an example 
in Figure [5] Within the range of nearest-neighbour coupling, each side of the outmost 3-fold gold atoms in 
the extended molecule is coupled to the gold electrode through nine gold atoms in the first layer and seven 
gold atoms in the second layer. All the gold-gold distance is fixed on 2.88A, which is the bond length of 
body gold. With the above procedure, we construct all the five gold-molecule-gold systems and show them 
in Figure |31 For simplicity, all the structures are treated symmetrically. 

3.2 The alignment of the Fermi level 

Another key point for calculating the conductance in molecular junctions is the position of the Fermi level Ep. 
As what Patrone has pointed out [16], UPS experiments can give direct information about Ep — Ehomo, 
which is an important energy reference in the molecular transport calculation. Because we calculate the 
electrodes and extended molecule separately, we have to face the difficulty that how to unify the zero 
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energy points of two different methods. Fortunately, introducing the gap of Ep — Ehomo informed by UPS 
data, we surpass this problem. With energy translation on the Papaconstantopoulpos diagonal Hamiltonian 
matrix elements [23], we can adjust the Fermi level in gold electrodes and locate it in the suitable place 
within molecular levels obtained by DFT calculations. Thus, the adjustment of the relative energy scale 
with respect to the experimental data of UPS [15, 17-19] becomes one of the meaningful procedures in our 
calculations. In Table 1 we show the calculation results of all the five different molecular junctions compared 
with the experimental values [6-8]. To our satisfaction, we found in our calculation that the conductance 
measured in STM experiments qualitatively coincide well with the UPS data via the bridge of Ep — Ehomo- 

As what we have shown in Table 1, after considering the UPS data of Ep — Ehomo 0.6eU [18], the 
calculation of BDT's transmission near the Fermi level (T(Ef) = 0.0086) agrees well with the experimental 
result 0.011. UPS experiments also show that the position of HOMO is larger than 5eV from the Fermi 
level in octandecanethiol monolayers on Au [15, 19]. For the similar electronic structures of three saturated 
alkancdithiol molecules (HDT, ODT, DDT), here we assume 5.5eU as an median value for the calculation of 
all the three molecules. The results of 0.85 x 10- 3 (HDT), 3.17 x l()- 4 (ODT), 2.62 x lO' 5 (DDT) could be 
qualitatively in accordance with the experimental results 1.2x 10" 3 (i7DT), 2.5x 10~ 4 (ODT), 2x 10" 5 (L>DT) 
[6-8]. Although we do not find direct UPS experimental results about BDMT, the measurement of 4,4- 
(ethynylphenyl)-l-bcnzenethiol [17], which has a comparable electronic structure with the incorporation 
of phenyl and carbon-chain group like BDMT, can give us useful enlightenments. Unlike the small value 
(0.5 ~ 0.6eU) in pure 7r-bonded molecules [16, 18] and the large value (> 5eV) in saturated carbon-chain 
molecules [15,19], the above mixed electronic structure has a 1.9eV gap of Ep —Ehomo [17]. It provides the 
prompt on estimating the Ep — Ehomo value of BDMT. Here we use 2.1eV and find that the calculation 
result 6.79 x 10~ 4 is in good agreement with the experimental data 6 x 10~ 4 . In addition, we consider 
the influence on the calculation results due to the uncertainty in UPS experiments. In Table 1, we change 
Ep — Ehomo informed by UPS data with a value of ±0.2eU. Our calculation shows that the results are not 
affected qualitatively. It indicates that introducing UPS experimental information might provide a relatively 
reliable energy reference for qualitative calculations on molecular conductance under equilibrium. 
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Table. 1 The used Energy gap of Ep — Ejjomo informed by UPS data [15, 17-19] and 
the comparison of calculated molecular conductances with the experimental results [6-8] 
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3.3 Analysis on density of states (DOS) 

The above results indicate the particular molecular electronic structures would associate with the alignment 
of the Fermi level in the gap between HOMO and LUMO. Basically, the energy position of the Fermi level 
and the degree of derealization in molecular junctions are considered as two major factors which impact the 
transport properties contributed from molecular orbitals [37] . In order to make a better comprehension on 
such an issue, we do some analysis on DOS of molecular junctions. For the consideration of charge transfer 
between gold electrodes and the molecule, here we use a central cluster (Au- molecule- Au) : one gold atom on 
each side of the molecule is included in our calculation. Mulliken's population analysis shows only a small 
charge transfer (0.12 ~ 0.17e) between the central cluster and the remained gold electrode parts in the above 
gold-molecule-gold models. So we believe that the main charge transfer has been properly described with a 
single gold atom connected to each side of the molecule. 

In Figure 0] we show the transmission under zero bias and the DOS of different parts in Au-BDT-Au 
as a function of energy. We find that the contribution from Au is significant in the gap between molecular 
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levels. In the junction of BDT, about 75% of DOS near the Fermi level is from gold. And in the other 
four junctions, the percentage of gold's contribution is around 90%. It means that gold electrodes might 
play an important role in molecular transport. When the coupling between molecules and electrodes is weak 
(in gold-HDT-gold, gold-ODT-gold, gold-DDT-gold junctions), the DOS of gold is mostly localized around 
the gold-molecule interface and leads to a tunneling behavior. Cui et.al. have discussed thoroughly on this 
subject of the transport in alkancdithiol molecular junctions by a tunneling model [38]. Their results are in 



accordance with our calculation shown in Figure 5(b) In addition, Kami's calculation suggests that transport 
accurs through HOMO [39] while our results show transport might accur through LUMO in alkanedithiol- 
gold junctions. On the other hand, the conductance in gold-BDT-gold junction is larger than the one in 
alkanedithiol junctions. The existence of n bond results in a strong hybridization of molecular states and 
gold's state, which implies a more complex transmission mechanism far beyond tunneling. 

We can also find additional information about transmission in the projected DOS (PDOS) on different 
molecular orbitals of Au- molecule- Au clusters. Here we reference the scheme of Heurich et. al. [37]. We 
choose the clusters including Au atoms for the decomposition on molecular orbitals. Such a scheme can 
directly reflect on the change of molecular orbitals induced by the adsorption on electrodes. We define new 
molecular orbitals after the hybridization of molecular states and gold' state in Au-molecule-Au clusters, 



which are meaningful but different from the ones in isolated molecules. Figure 6(a) Figure |7(a)| and Figure 



8(a) respectively show the PDOS of several representative molecular orbitals in Au-BDT-Au, Au-BDMT-Au 



and Au-HDT-Au. The corresponding orbital-shape plots are shown in Figure 6(b) Figure 7(b) and Figure 



8(b) According to these figures, the contribution of a special molecular orbital (MO) to DOS strongly 



depends on its charge-density distribution on this MO. The extended and well-coupled MOs, such as the 
HOMO, LUMO and LUMO+1 in Au-BDT-Au, have a significant weight in DOS. After coupled with gold 
through 6s atomic orbital, these molecular orbitals extended widely (> 6eV). The HOMO and LUMO+1 
are even partly hybridized. The derealization of 7r-bonded electronic structure is the main reason for the 
existence of such MOs. On the contrary, the localized PDOS always derives from the confined charge-density 



on this MO. The orbital shape of Au-HDT-Au in Figure 8(b) shows a distinct feature of localization about a- 
bonded molecules. Figure0also shows an obvious influence of -(CH2)- structure on transmission. It isolates 
the coupling between electrodes and the 7r bond in BDMT and results in a localized DOS distribution of 
most MOs. 



4 Conclusion 

In this work, the conductances of a series of gold-dithiol molecule-gold junctions are calculated by matrix 
Green's function method combined with hybrid tight-binding/ density functional theory. The atomically- 
contacted extended molecule model is used for the calculation. After determining the position of the Fermi 
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level by the energy reference according to the results from UPS experiments, the final calculation of molecular 
conductances qualitatively agree with the experimental values obtained by Tao and co-workers [6-8] . In our 
calculation, the influence by the uncertainty of UPS data is also considered. Our results show that the first- 
principles method based on Green's function method can qualitatively predict the conductance properties 
near low bias under the instruction of UPS experiments. 
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Figure 1: The optimized geometry for Au3-BDMT-Au3. The optimization method is B3LYP/6-311G** + LANL2DZ. 



Figure 2: The example of extended molecule: Au4-BDMT-Au4 




Figure 3: Au20-molecule-Au20 clusters used in the calculations of conductances with B3LYP/6-311G** + LANL1MB. 
The molecules are benzenedithiol (BDT), benzenedimethanethiol (BDMT), hexanedithiol (HDT), octanedithiol (ODT) and 
decanedithiol (DDT). 
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Figure 4: (a)The transmission coefficient (b) The contribution to total density of state (TDOS) of Au-BDT-Au from different 
parts. The blue solid line is the part from phenyl ring. The red dashed line is the part from sulfur. The magenta thick line is 
the part from gold. In the inset of the amplified figure of DOS, we can see the DOS near the Fermi level. 
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Figure 5: (a) Transmission coefficient of Au-BDMT-Au junction, (b) Transmission coefficient of Au-HDT-Au, Au-ODT-Au, 
Au-DDT-Au junctions. 
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Fi gure 6i (a) The projected density of states (solid line) of six molecular orbitals of Au-BDT-Au. The dashed line is the 
total density of state of Au-BDT-Au. (b) The orbital shapes of the corresponding molecular orbitals in (a) 
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Fi gUT6 7; (a) The projected density of states (solid line) of six molecular orbitals of Au-BDMT-Au. The dashed line is the 
total density of state of Au-BDMT-Au. (b) The orbital shapes of the corresponding molecular orbitals in (a) 
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Fi gure 8: (a) The projected density of states (solid line) of six molecular orbitals of Au-HDT-Au. The dashed line is the 
total density of state of Au-HDT-Au. (b) The orbital shapes of the corresponding molecular orbitals in (a) 
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